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ABSTRACT 

We present 2 — 5 @m adaptive optics (AO) imaging and polarimetry of the famous hypergiant stars 
IRC +10420 and VY Canis Majoris. The imaging polarimetry of IRC +10420 with MMT-Pol at 2.2 /im 
resolves nebular emission with intrinsic polarization of 30%, with a high surface brightness indicating 
optically thick scattering. The relatively uniform distribution of this polarized emission both radially 
and azimuthally around the star confirms previous studies that place the scattering dust largely in 
the plane of the sky. Using constraints on scattered light consistent with the polarimetry at 2.2 /im, 
extrapolation to wavelengths in the 3 — 5 /im band predicts a scattered light component significantly 
below the nebular flux that is observed in our LBT/LMIRCam 3 — 5 /im AO imaging. Under the 
assumption this excess emission is thermal, we find a color temperature of ~ 500 K is required, well in 
excess of the emissivity-modified equilibrium temperature for typical astrophysical dust. The nebular 
features of VY CMa are found to be highly polarized (up to 60%) at 1.3 /im, again with optically 
thick scattering required to reproduce the observed surface brightness. This star’s peculiar nebular 
feature dubbed the “Southwest Clump” is clearly detected in the 3.1 /im polarimetry as well, which, 
unlike IRC+10420, is consistent with scattered light alone. The high intrinsic polarizations of both 
hypergiants’ nebulae are compatible with optically thick scattering for typical dust around evolved 
dusty stars, where the depolarizing effect of multiple scatters is mitigated by the grains’ low albedos. 
Subject headings: stars: supergiants - stars: individual (IRC +10420, VY Canis Majoris) - stars: 

circumstellar matter - stars: winds, outflows - instrumentation: polarimeters - 
techniques: high angular resolution 


1. INTRODUCTION 


Hypergiant stars are located near the upper limit of 
luminosity in the HR Diagram, with typical luminosi¬ 
ties of L* ~ 5 x 10 5 Lq and effective photosphere radii 
of several AUs. They exhibit very h igh mass loss rates 
of as much as 10 -4 M 0 yr _1 (Danchi et al. 1994) with 


some experiencing even greater discrete eruptions. Here 
we present new adaptive optics (AO) infrared observa¬ 
tions of two famous hypergiant stars IRC +10420 and 
VY Canis Majoris. 

IRC +10420 is a yellow hypergiant, in the rare short¬ 
lived evolutionary stage experienced by massive stars 
transiting to or from the red supergiant stage. Hub¬ 
ble Space Telescope (HST) visual imaging revealed a 
very complex circumstellar nebula with numerous con- 


HST long-slit spectroscopy 
2002) combined with second epoch 
HS'l' imaging found”that we view the sta r nearly pole- 


(Humphreys et al. 

1997 

(Humphreys et al. 

2002 


on,_k>oking down on its equatorial plane (Tiffany et al. 
2010). Recent VLTI observations of neutral and ionized 


gas arou nd IRC +104 20 have confirmed this geometry 
(Oudmaijer & de Wit 20"l3|). 

VY ClVla is a cool, red hypergiant. Multi-epoch HST 
imaging and polarimetry found extensive episodic mass 
loss with no preferred axis of symmetry (Humphreys 


Observations reported here were obtained at the MMT Ob- 

servator^ajoint^a+lh^omhe^Smithsonmjmnstitutior^ncmhe 

University of Arizona. 


et al. 2007 Jones et al. 2007). Ground-based 2 - 5 
pt m AO imaging with LlVllRCam ( Skrutskie et al.|[~2010 ) 
on the La rge Binocular Telesc ope (LBT|/] previously re¬ 


ported by Shenoy et al. (2013) found VYuMa’s peculiar 
“Southwest Clump” to be a particularly dense, optically 
thick condensation of the ejecta. Most recently, ALMA 
sub-millimeter observations have found further discrete, 
dense ejecta in the very close environment of the star 
(|Richards et al.|2014| |0’Gorman et al.|2015|) . The mech¬ 
anisms lor such discrete mass loss are not understood. 
Magnetic field activity an alogous to coronal mass ejec¬ 
tions may be responsible (Humphreys et al. 2007). Al¬ 
though recent XMM — A lewton X-ray observations of 
VY CMa placed constrains on the magnetic field strength 
at the star’s surface, the star may have simply been in 
state of lower m agnetic activit y than during previous 
mass loss events (Montez et al. 2015). 

Because of the high contrast ratio between the star’s 
light profile and the surface brightness of nebular mate¬ 
rial, separating infrared nebular emission from the star 
presents considerable difficulty. Imaging polarimetry is 
a powerful technique for studying nebular emission in 


2 The LBT is an international collaboration among institutions 
in the United States, Italy and Germany. LBT Corporation part¬ 
ners are: The University of Arizona on behalf of the Arizona uni¬ 
versity system; Istituto Nazionale di Astrofisica, Italy; LBT Beteili- 
gungsgesellschaft, Germany, representing the Max-Planck Society, 
the Astrophysical Institute Potsdam, and Heidelberg University; 
The Ohio State University, and The Research Corporation, on be¬ 
half of The University of Notre Dame, University of Minnesota and 
University of Virginia. 
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TABLE 1 

Summary of Observations 


Date 

Instrument 

Filter(s) 

On-Source Time 

(UT) 


(/im) 

(<0 


IRC 

+10420 


2011 May 25 

LMIRCam 

PAH1 (3.29) 

186 

1! 

11 

PAH2 (3.40) 

220 

1! 

11 

L' (3.83) 

154 

1! 

11 

M (4.9) 

K' (2.2) 

107 

2012 Sep 25 

MMT-Pol 

254 


VY CMa 


2013 Oct 22 

MMT-Pol 

J' (1.3) 

36 

2014 Jan 16 

! ! 

3.1 

320 


this regime. Dust in the nebula scatters light from the 
star and polarizes it. If the star’s light is unpolarized or 
weakly polarized compared to the fractional polarization 
of the scattered star light, then imaging polarimetry can 
cleanly separate the faint polarized light from the star’s 
dominating light. We present new observations of these 
two hypergiant stars made with MMT-Pol, the 1 — 5 /im 
imaging polarimeter custom d esigned for the 6.5 m M MT 
observatory on Mt. Hopkins (Packham et al. 2012) and 
with LBT/LMIRCam on Mt. Graham, AZ . The obser- 
vations are summarized in Table 1. 


2. OBSERVATIONS AND RESULTS 

2.1. IRC + 10420 : 2.2 i±m AO Imaging Polarimetry 

We obse rved IRC +10420 on 2012 Sep 26 (UT) with 
MMT-Pol (Packham et al.||2012|) during its commission¬ 
ing run. MMT-Pol operates at the Cassegrain focus, 
where the AO-corrected beam from the secondary mir¬ 
ror enters the instrument with no prior off-axis reflections 
that might introduce spurious polarization. A rectangu¬ 
lar mask in the focal plane provides a 20" x 40" field 
of view. The beam passes through a half-wave retarder 
(or half-wave plate, hereafter HWP), after which a Wol¬ 
laston prism splits the light into orthogonal components. 
These two components are referred to as the ordinary 
and extraordinary rays, or o and e rays. Rotating the 
HWP from 0° to 45° swaps the o and e rays. Compar¬ 
ison of the two rays allows Stokes parameter Q to be 
extracted while canceling out any differences in through¬ 
put between the two rays and/or differences in detector 
response where the rays are detected. Similarly, rotating 
the HWP from 22.5° to 67.5° yields Stokes parameter U. 
MMT-Pol’s pixel scale is 0.043" pix -1 . 

We imaged IRC +10420 through the K 7 (Ao = 2.20 /im, 
A A = 0.12 /im) filter at two dithered positions separated 
by 20" on the detector array. At each dither position im¬ 
ages were taken as correlated double-sample (CDS) pairs, 
with the subtracted difference of each CDS pair yielding 
a single 1.33 s exposure. For each HWP position an¬ 
gle 50 exposures were taken. For image quality and flux 
calibration, a polarized standard star (HD 29333) and 
unpolarized standard star (HD 224467) were observed. 
The FWHM of the PSF was 0.2". The flux calibration 
factors (ADUs to W cm -2 ) computed using these two 
standards agree to within 5%. The instrumental polar¬ 
ization determined from the unpolarized standard star 


was 0.47% ± 0.07%. 

Each 1.33 s exposure was examined individually. On 
occasion when an increase in the width of the star pro¬ 
file indicated potential loss of AO lock, those exposures 
were discarded. The remaining useable frames were mean 
combined to yield a single image in each HWP position 
for each of the two dithers. The sky background and ar¬ 
ray structure were removed by subtracting dithered im¬ 
ages; this subtraction removed many hot pixels as well. 
Remaining hot pixels were replaced with the median of 
the surrounding 8 pixels. The o-ray and e-ray images 
were extracted from these images. The central point 
source of IRC +10420 in each exposure saturates the 
images out to a radius of < 0.4". For each dither po¬ 
sition, the 8 images (2 per HWP position) were aligned 
using the IRAFq] XREGISTER task. The images were 
then smoothed with a Gaussian of FWHM = 0.2" (the 
beam size as estimated from the PSF star). 

The normalized Stokes parameters q = Q/I and u = 
U/I were computed using the ratio method of [Tinbergen 
(1996): 


Q = 


R q - 1 

Rq + l 


, for R q = 



Ru- 1 
Ru +1 


, for R u = 




22 


ih/i't 


67 


(i) 


( 2 ) 


where the o and e superscripts on each of the intensities 
/ refer to the ordinary and extraordinary rays emerging 
from the Wollaston prism, and the subscripts refer to 
the four HWP orientation angles 00.0°, 45.0°, 22.5°, and 
67.5°. Dividing the intensities to form ratios R q and R u 
minimizes potential differences in the Wollaston prism’s 
relative throughput on the two rays, relative differences 
in the response of the array locations where the two rays 
are detected and variations of sky transmission during 
the time between when the images are taken through the 
paired orientations of the HWP. The normalized Stokes 
parameters q and u computed with this ratio method are 
algebraically equivalent to the standard definitions: 

_ Q _ 1 0 — I90 _ U _ I45 — Cl 35 , x 

q ~ I ~ Jo + J90 5 U ~ I ~ h 5 + 1 135 1 ’ 

where are the intensities transmitted through a linear 
polarizer (analyzer) with its transmission axis at angles 
(j) = 0°, 90°, 45° and 135°, and / is the total intensity. 
The fractional polarization p and position angle (PA) 9 
are then computed with: 


p= \/q 2 +u 2 ; 6 = ^ tan 1 (4) 

The polarized intensity is the product of the fractional 
polarization and total intensity: Ip = p • I. An S/N cut 
of 4-cr or 8-cr with respect to background fluctuations in 
polarized intensity off the source was selected as the cut¬ 
off for the display of polarization vectors. A 4-sigma cut 
is equivalent to de-biasing by a factor of less than 0.03 


3 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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(see Equation (A3) of Wardle & Kronberg (1974])), which 
has no impact on our science. 

The position angle 0 computed using Equation (4) is in 
array coordinates. To convert to degrees East-of-North 
on the sky requires the addition of a calibration offset an¬ 
gle A#, which is normally determined from observations 
of one or more polarized standard stars. Observations of 
the polarized standard star HD 29333 during the com¬ 
missioning run did not provide a meaningful calibration 
offset due to non-photometric conditions at the time the 
standard was observed and inadequate integration time. 
For our observations of IRC +10420 the uncorrected po¬ 
larization vectors’ orientation is physically implausible. 
Dust distributed around and illuminated by a single cen¬ 
tral source is expected to show a centro-symmetric scat¬ 
ter ing^attem, and previo us observations of IRC +10420 
by Kastner & Weintraub (1995) find such a pattern. An 
offset of A Ok' = 58 u when added to each of our uncor¬ 
rected polarization vectors achieves a centro-symmetric 
pattern in all directions around the star. We have there¬ 
fore adopted this offset value. 

In Figure [lja) we display our resulting p raw , the raw 
polarization around IRC +10420 prior to separation into 
stellar and nebular components. By “raw” we mean un¬ 
corrected for any underlying flux from the star in the 
wings of the PSF. With the adopted PA offset, the raw 
polarization is consistent with the centro-symmetric pat¬ 
tern characteristic of scattering of light from the central 
star by circumstellar dust. The raw percentage of po¬ 
larization increases with radius, rising to a maximum of 
15% at 1.7". This suggests the central star profile dilutes 
a higher intrinsic polarization of the nebula; this effect 
is corrected for in Section §3.1 below. In Figure [TJb) 
we display a flux-calibrated image of the raw intensity 
I raw with the color scale in mag arcsec -2 for compari¬ 
son with the contours in 1(a). In Fig. [ljc) we display 
the product of these two images, which is the polarized 
intensity Ip LW = p raw • I raw . While the raw intensity 
is dominated by the star and thus drops off smoothly 
with radius, the raw polarized intensity image shows a 
region of relatively flat, extended emission from 0.5" out 
to a radius of ~ 2.5". This corresponds to emission from 
extended nebular material, not the stellar PSF. 

2 .2. IRC + 10420 : 3-5 jinn AO Imaging 

IRC +10420 was observed on 2011 May 25 UT dur¬ 
ing the commissioning of LMIRCam, the 2 — 5 /im 
high-res olution camera on the L arge Binocular Telescope 
(LBT) (Skrutskie et al. 2010). During commissioning 
a single 8.4 m primary mirror was used, in conjunction 
with the deformable AO secondary. This system achieves 
near-diffraction limited imaging and high sensitivity in 
the thermal infrared due to minimizing the number of 
reflecting surfaces and keeping all mirrors after the ter¬ 
tiary mirror at cryogenic temperatures. IRC +10420 was 
observed through two narrow-band filters: PAH 1 (Ao = 
3.29 /im, AA = 0.02 gm), PAH 2 (A 0 = 3.40 /im, AA = 
0.02 /im) as well as two broad-band filters: L' (Ao = 3.8 
/im, AA = 0.3 /im) and M (Ao = 4.9 /im, AA = 0.3 /im). 
LMIRCam has a field of view approximately 10" x 10", 
with a pixel scale of 0.011" pixel -1 . The images of IRC 
+ 10420 were dithered for background sky subtraction. 
For the two PAH filters, the star 77 Aquilae provid es a 
PSF (77 Aql is point-like in the 3-8 /im range (Barmby 


et al. 2011|)), while for the L' filter the star BD +35 2435 
is the PSF. 

In Figure |2] we display the four LMIRCam images of 
IRC +10420 along with the PSF stars for comparison. In 
the bottom row of the figure, we have subtracted a scaled 
Moffat profile fitted to the PSF at each wavelength. For 
the M filter we radially scaled the L' PSF profile by the 
ratio of A m to A//. In all four filters IRC +10420 is 
substantially extended compared to the PSF. 

We used two methods to determine a mean flux cal¬ 
ibration factor in each filter. For the first method we 
equated the integrated intensity in ADUs per second out 
to the 3- <7 level in each image (i.e., entire source = star + 
extended emission) with the specific flux F v from the ISO 
SWS spectrum for IRC +10420. For the second method, 
we calibrated the ADUs per second against those of the 
observed standard stars. For the flux of the standard star 
77 Aql at the PAH 1 & PAH 2 filter wavelengths we inter¬ 
polated between its reported magnitu des in th e 2 M ASS 
Ks filter and IRAC Band 1 (Marengo et al.||2010). For 
the flux of the standard star BD +35 2435 at 17 we inter¬ 
polated between its WISE W1 and W 2 magnitudes. In 
the absence of a standard star observation in the M fil¬ 
ter, we repeated the first method with comparison to IRC 
+ 10420’s magnitu de at M from the large, single aperture 
photometry of Jones et al. (1993). For each filter we use 
the mean calibratlonJactors from the two methods, with 
estimated uncertainties of 15%, 10%, 10% and 5% in the 
PAH 1 , PAH 2 , L' and M filters respectively. 

2.3. VY CMa: J' (1.3 /am) and 3.1 /im Imaging 
Polarimetry 

We observed VY CMa with MMT-Pol on 2013 Oct 22 
UT through a narrow-band J' (Ao = 1.3 /im, AA = 0.1 
/im). Six sets of dithered images were taken with a 1 
s exposure time per HWP position. Data reduction was 
performed as discussed above for the MMT-Pol images of 
IRC +10420. The unpolarized standard star HD 224467 
provided the PSF and flux-calibration. The offset A 6j> 
added to all polarization vectors was chosen to yield a 
centro-symmetric scattering patter n consisten t with the 
HST visual polarimetry in Jones et al. (2007). We esti¬ 
mate our photometric uncertaintyDn the nux-calibration 
from the two standards at 7%. We display our J' po¬ 
larimetry in Figure [3j and a comparison to the HST 
visual polarimetry in Figure |4j 

We subsequently observed VY CMa through MMT- 
Pol’s narrow-band 3.1 /im filter (AA = 0.1 /im) on 2014 
Jan 16 UT. The PSF and flux-calibration were obtained 
from observations of HD 104624 and HD 29333. We esti¬ 
mate our photometric uncertainty in the flux-calibration 
from these two standards at 5%. For the 3-5 /im range, 
MMT-Pol uses a retarder which provides a retardance 
of 180° (half-wave) at 3.6 /im. For 3.1 /im light, the 
waveplate’s retardance is 209°. Defining efficiency as the 
ratio of the input polarized intensity to the measured po¬ 
larized intensity, the difference in retardance lowers the 
efficiency of the polarization measurement at 3.1 /im by 
a factor of 0.94. We have corrected the raw polarized in¬ 
tensity by this 6 % factor during the reduction of the 3.1 
/im images. The data reduction is otherwise the same as 
for IRC +10420, except we applied a S/N cut-off of 8 -cr 
in polarized intensity to mitigate faint artifacts from the 
heavily over-exposed central star. A polarization vector 
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IRC+10420: K' (2.2 pm) 


p raw (raw polarization) 



A RAH 


r aw (raw total intensity) 



4 3 2 1 0 -1 -2 -3 -4 

A RA (") 


Ip aw (raw polarized intensity) 



A RA (") 



(a) (b) (c) 

Fig. 1. — MMT-Pol observations of IRC +10420 at A = K' (2.2 pm). North is up and East is left. The center is masked out to a radius 
of 0.4" where the images saturated. The circle in the lower left is the beam size (0.2" = FWHM of PSF). (a): The vectors are the raw 
polarization p raw for a cut in polarized intensity of 4-a with respect to background fluctuations off the source. The overlaid contours show 
the raw total intensity (total = unpolarized + polarized) from 9 to 6 mag arcsec -2 , progressing inwards in steps of —1 mag arcsec -2 . A 
portion of the image was masked in the south where faint artifacts which appear to have been caused by internal reflections of the heavily 
overexposed central point source within the instrument optics cause spurious polarization. The vectors show a centro-symmetric pattern 
characteristic of scattering by circumstellar dust. The raw polarization rises to about 15% at a radius of ~ 1.7" around most of the star. A 
length scale for the polarization vectors is given in the lower right of the image, (b): Raw total intensity I raw in mag arcsec -2 , with the 
same overlaid contours as left image, (c): Raw polarized intensity I™ w = p raw . jraw - n ma g arcse c -2 . The dashed-line annulus indicates 
the region of overlap-of-coverage (OCR) with the 3 — 5 pm LMIRCam observations (Fig. |5|. A colorbar for the surface brightnesses is on 
the far right. 
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Fig. 2. — IRC +10420: 3 - 5 pm adaptive optics images made 
with LBT / LMIRCam. Top Row: IRC +10420, with the 
filter for each column of images as indicated. Each FOV is 5" X 
5", with square-root scaling to emphasize the extended nebular 
emission. Middle Row: PSF stars. For filters PAH1 and PAH2 
the PSF is the star 77 Aql (FWHM = 0.10" in each filter), while at 
V the PSF is BD +35 2435 (FWHM = 0.11"). There is no PSF 
image for the M filter from that night. Bottom Row: Same as 
the top-row images, after subtracting a scaled Moffat fit to each 
filter’s PSF. To simulate a PSF for the M filter, the FWHM of the 
Moffat function used at L' was scaled by the ratio AThe 
top and middle row images are each stretched to the maximum 
pixel value at the center. The subtracted images of IRC +10420 
in the bottom row are stretched to the same maximum value as 
the top row images. The dashed-line annulus on the subtracted 
images is the region of overlap-of-coverage (OCR) with the 2.2 p m 
polarimetry. It spans from 0.5" (where the 2.2 pm polarimetry is 
outside the saturated stellar image) out to 1.0", the 3-cr level in 
the PAH1 (3.29 pm) and M images. 


offset has been added which achieves a centro-symmetric 
pattern. We display our result in Figure [5] 

3. ANALYSIS & DISCUSSION 

3.1. IRC +10420: Nebula’s Intrinsic Polarization at K' 

In_their lar ge-ap erture polarimetry of IRC +10420 
Jones et al. (1993) observed a fractional polarization 
of 17o at K' for the combined star plus any unresolved 
nebula. The sub-arcsecond resolution now possible with 
MMT-Pol allows us to resolve nebular emission for radii 
greater than about 0.4" (2000 AU), limited mostly by 
the exceedingly bright central star, which necessarily sat¬ 
urates in the integration times required to bring up the 
nebulosity. In our image the raw fractional polarization 
p raw varies with radius, rising with increasing distance 
from the star up to a peak of about 15% at 1.7" and 
then dropping off. This behavior indicates that flux in 
the PSF from the weakly polarized central star is dilut¬ 
ing a higher intrinsic polarization of the nebula nearer 
the star. 

To assess the nebula’s intrinsic polarization, we first 
estimate the nebula’s total (unpolarized + polarized) in¬ 
tensity I neb by subtracting a representative star profile 
from the raw intensity I raw . On the left-hand side of 
Figure [6] we plot the azimuthal average profile of I raw 
for the regions to the East and West of the star, ex¬ 
cluding the regions to the North and South where faint 
artifacts from the star are present. The star’s light pro¬ 
file /* is made using a Moffat function fitted to the PSF 
so that 7* can be extended out into the nebula. The fit¬ 
ted star profile has been scaled vertically so that it does 
not over-subtract by a radius of about 0.5". This scaling 
consistent with the difference in observed K magnitudes 
of the PSF star and IRC +10420, assuming K = 3.5 for 
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p raw (raw polarization) 



VY CMa: )' ( 1.3 /xm) 
f aw (raw total intensity) 



A RAP) 


Ip aw (raw polarized intensity) 
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(a) (b) (c) 

Fig. 3.— MMT-Pol observations of VY CMa at A = J' (1.3 /xm). North is up and East is left. The center is masked out to a radius of 0.6" 
where the star saturates, (a) : The vectors are the raw polarization p raw for a cut in polarized intensity of 4-a with respect to background 
fluctuations off the source. The overlaid contours show the raw total intensity (total = unpolarized + polarized). The outermost contour 
is 10 mag arcsec -2 , with contours progressing inwards in steps of —1 mag arcsec -2 . The extended shape of the outermost contour towards 
the North is due to a faint artifact which appears to have been caused by internal reflections of the heavily overexposed central point source 
within the instrument optics. In contra st th e extensio n to the South is real, coinciding with the location of the distinct component of the 
ejecta previously identified as Arc 2 by [Humphreys et al. (2007) (see Fig. [4ja) below). The polarization vectors show a centro-symmetric 
pattern characteristic of scattering by circumstellar dust, (b): Raw total intensity I raw in magnitudes arcsec -2 , with the same overlaid 
contours as left image, (c): Raw polarized intensity I™ w = p raw . jraw - m ma g arcsec -2 . A colorbar for the surface brightnesses is on the 
far right. 



ARA (") ARA (") ARA (") 


(a) (b) (c) 

Fig. 4.— VY CMa: (a) HST F1042M (1 /xm) image reproduced from Humphreys et al. I (120071) identi fying the NW Arc, Arc 2, S Knot, 
S Arc and SW Clump, (b) HST visual (0.66 /xm) polarimetry replotted using data trom ]Jones' et a l. ] 42007 ). (c) MMT-Pol J' (1.3 H 
polarimetry (this work). 


the latter (Oudmaijer et al.||2009). The difference of the 
observed and star profiles is deemed to be the nebula’s 
total intensity: I neb = I raw — /*, which is depicted in 
Figure [6] with open circles. 

For the nebula’s polarized intensity Ip eb we assume 
that all of the raw polarized emission in Figure [ljc) is 
entirely from the nebula: Ip <C Ip eb so that Ip LW —» 
Ip eb . This only slightly overestimates Ip eb within a ra¬ 
dius of about 1" and becomes increasingly accurate as 
radius increases. For example with a 1% polarized star 
whose profile follows the PSF profile, by a radius of 1" 
the star’s polarized intensity Ip is ~ 10% of Ip aw and 
falls off rapidly as the star’s light profile descends going 
farther into the nebula region. The radial profile of Ip eb 
is plotted with open squares in the upper right of Figure 
[6] along with the nebula’s total intensity I neb . The ratio 
of these intensities is used to compute the nebula’s in¬ 
trinsic fractional polarization, which by definition is p neb 


= Ip eb / I neb . This intrinsic polarization p neb is plotted 
with solid black squares in the lower right of Figure 6j 
This removal of the star’s light profile indicates that the 
nebula’s intrinsic polarization is about 30% over a broad 
radial range from 0.5" to ^ 2 ". 

Based on radial vel ocity and proper motion measure¬ 
ments, Tiffany et al. (2010) found that the ejecta seen 
in HST visual images is moving mostly in the plane of 
the sky, and concluded that IRC+10420 is nearly pole-on 
and we are looking down onto an equatorial dust distri¬ 
bution. The relatively high fractional polarization (30%) 
at K' that we find for the nebulosity surrounding the star 
is consistent with this model. For pure Rayleigh singly 
scattered light (size parameter 27ra/A <C 1), polariza¬ 
tion as high as 100% could in principle be expected for 
optically thin dust located exactly in the plane of the 
sky with respect to the star (scattering angle 0 = 90°). 
Realistically the maximum possible polarization is likely 
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p raw (raw polarization) 



VYCMa: 3.1 /im 


F aw (raw total intensity) 
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Ip aw (raw polarized intensity) 
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Fig. 5.— MMT-Pol observations of VY CMa at A = 3.1 p.m. North is up and East is left. The center is masked out to a radius of 1.0 // 
where the star saturates, (a) : The vectors are the raw polarization p raw for a cut in polarized intensity of 8-cr with respect to background 
fluctuations off the source. The overlaid contours show the raw total intensity (total = unpolarized + polarized). The outermost contour 
is 5 mag arcsec -2 , with contours progressing inwards in steps of —1 mag arcsec -2 . (b): Raw total intensity I raw in magnitudes arcsec -2 , 
with the same overlaid contours as left image, (c): Raw polarized intensity I™ w = p raw . jraw - in ma g arcsec -2 . A colorbar for the 
surface brightnesses is on the far right. The raw polarization at the location of the SW Clump is 20%. 


lower for one or more reasons. For a distribution of grain 
sizes that includes larger grains with 2i:a/\ on the order 
of 1, Mie theory applied to spherical grains indicates the 
maximum polarization decreases below 100%. Consider 
for example a Gaussian distribution of grain sizes with 
mean radius a = 0.3 fi m and standard deviation a a = 
0.15 fi m, with astronomical silicat e optical constants a t 
2.2 fi m of n = 1.72, k = 0.035 (Draine & Lee 1984) 
Computing the maximum polarization for such a distri 


bution with the BHMIE subroutine (Bohren & Huffman 
1983) yields a maximum polarization of ^ 60% occur 


ring at a scattering angle slightly behind the plane of 
the sky. Alternately, for a distribution of scattering an¬ 
gles about 90° (a flared disk, for example) and multiple 
scatters, one would expect a lower fractional polarization 
than the ideal maximum. Since the nebulosity as seen in 
polarized intensity is distributed in a nearly circular pat¬ 
tern around the star and has a near constant fractional 
polarization, it can not be in a disk or torus that is tilted 
to the line of sight. For such a disk or torus the projected 
shape would be elliptical and the scattering angles would 
vary with azimuthal angle. The only other possible ge¬ 
ometry that could explain the observed morphology and 
polarization of the nebulosity surrounding IRC-f 10420 
is a thin cone opening towards (or away) from us at a 
uniform scattering angle and surface brightness with dis¬ 
tance from the star. We find this to be implausible in 
light of the radial velocity measurements, although not 
ruled out by our polarimetry. 


3.2. VY CMa: Intrinsic Polarization of Nebular 
Features at J (1.3 /am) and 3.1 pm 

The polarized intensity images of VY CMa in Figures 
|3jc) & |5|c) reveal several distinct featu res of its complex 
nebu la!~Using the names designated by [Humphreys et ah] 
(|2007|) in their HST F1042M (1 /im) image (reproduced 
in Figure [4a]), at J' we observe the Northwest Arc, Arc 
2, the Southwest Clump, the South Knot, and the South 
Arc. The first two features are located sufficiently distant 
from the star that their raw polarization is taken to be 


their intrinsic polarization. The raw polarimetry at J' 
(1.3 /im) shows the NW Arc is ^35% polarized, while 
Arc 2’s polarization rises as high as 45%. These and 
subsequently discussed values are summarized in Table 
2 . 


For the SW Clump, S Knot and S Arc, the intrinsic po¬ 
larization of each is diluted by the star to varying extents. 
Applying the same procedure as in §3.1, we assume the 

raw polarized intensity is solely from the nebular feature 
b 


so that I™p 


I\ e p- We estimate the feature’s total 


intensity I^ eb by subtracting from the raw intensity the 
profile of the star’s light, which we represent using a pho¬ 
tometric cut at position angle +148° East of North. Each 
nebular feature’s intrinsic polarization is the ratio: p^ eb 
= / I% eb . For the SW Clump as shown in Figure 

[7^a) for example, this raises its J' raw fractional polar¬ 
ization slightly to an intrinsic polarization of 40%. For 
the S Knot and S Arc the same procedure raises the raw 
polarization of each from 40% to an intrinsic polarization 
of as much as 60%. 

At 3.1 /im the SW Clump is the sole feature observed 
in polarized intensity (Figure [5 'c)). Per Figure [7^b) 
subtracting the star profile in that filter indicates the 
SW Clump’s intrinsic polarization is a factor of about 
2 higher than the raw polarization, increasing from 20% 
to ~ 40%. The SW Clump is a particularly curious fea- 
ture in the ejecta of VY CMa. Humphreys et al. (2007) 
found it to be moving slowly away from the star with 
its motion largely in the plane of the sky. Previously 
reported AO imaging with LMIRCam has demonstrated 
that the Clump’s 2 — 5 /im emission must be due largely 
to diffusely reflected (i.e., optically thick) scattered light 
(Shen oy et al.||2013 ). 

A sili cate gram model fitt ed to the Clump’s 2 — 5 /im 
flux by Shenoy et al. (2013) indicated a dust mass lower 
limit of 5 x 10 -5 M 0 , with a temperature constrained to 
be between 80 — 210 K. Interestingly, subsequent ALMA 
mi llimeter-range (mm) obs ervations of VY CMa reported 
by O’Gorman et al. (2015) did not detect the Clump at 
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Fig. 6 .— Estimating the intrinsic polarization of IRC +10420’s nebula at K 7 (2.2 /im). Left: The filled circles are the azimuthal average 
of the raw intensity I raw to the East and West of IRC +10420, using radial bins equal to the beam size (0.2 // ). The dashed line is a 
Moffat profile fitted to the PSF, scaled to represent the star’s profile /*. The difference (open circles) is deemed to be the nebula’s average 
profile I neb . Upper Right: A zoomed view of the nebula’s average profile I neb (open circles). The open squares are the average profile 
of the raw polarized intensity I™ w (see Figure ^c)), which is then assumed to be solely nebular emission: —> I™ eb . Lower Right: 

The nebula-only polarized intensity is then divided by the calculated nebula’s total intensity to yield the nebula’s intrinsic polarization 
p neb = I™ eb j I neb (filled squares), compared to the raw polarization p raw (triangles). This removal of the star’s light profile shows that 
on average the nebula’s intrinsic polarization is at least a factor of 2 higher than the raw polarization and relatively constant from 0.5 // — 2". 




(a) (b) 

Fig. 7. — VY CMa: Photometric cuts through the SW Clump at J' (1.3 jtim) (left column) and 3.1 /im (right column). The SW Clump’s 
proximity to the star results in the star’s light profile diluting the Clump’s intrinsic polarization. Top row of (a) &; (b): For each 
indicated filter the thick solid line is the raw (total) intensity along a cut through the Clump at a position angle of —135° E of N. The thin 
solid line is a cut at a position angle of +148° E of N, which is taken to represent the star’s light profile into the region of the Clump. The 
dashed line is the difference of the two. Bottom row of (a) &; (b): For each filter the thick solid line is the raw fractional polarization 
p raw along the same cut through the location of the Clump. The Clump’s intrinsic fractional polarization p cl is estimated from scaling 
up p raw by the ratio of I raw j I cl from the plot above. 
























































321 & 658 GHz. For the lower-limit dust mass computed 
from the 2 — 5 fi m scattered light, those authors find 
it would require the SW Clump to have an unusually 
low temperature (<20 K) to be undetected by ALMA. 
Alternately, given the 3 -cr sensitivity limit of their obser¬ 
vations, for a temperature between 80 - 210 K they find 
a dust-mass upper limit of 3 — 7 x 10 -6 M 0 within a 
single 0.229" x 0.129" elliptical beam. In the infrared 
the Clump spans an area roughly ^ 1 " x ~ 0.5". For 
the area spanned by the beam size of the mm observa¬ 
tions, the mass lower limit from scattered light would 
decrease by a factor of about 3. Although that brings 
the infrared and mm mass estimates to within about a 
factor of 2 , some discrepancy between these mass esti¬ 
mates remains. One possible explanation is the assumed 
wavelength dependence of the dust emissivity. Typically 
the emissivity is parameterized with a power-law form 
where the emissivity = (A/Ao) - ^. For the^ diffuse 
ISM, 8 takes on the value of 1.6 (|Planck Collaboration 


et al. 2014). O’Gorman et al.| ( 2015| ) use a lower value 


of [3 = 0.7 based on their own spectral index fits and 
previous work on mm e mission from VY CM a and simi¬ 
lar la t e typ e stars (e.g. Knapp et al. (1993)). |Kamiriski 


et al. (2013) find (3 ~ — 0.5 for VY CM a/s 2/ 9 — 355 GHz 


continuum emission. They do not see the SW Clump in 
their continuum map, however they do detect it in their 
line maps of H 2 S (300.5 GHz) and CS (293.9 GHz). The 
absence of mm thermal dust continuum emission from 
the SW Clump could be due to a value of (3 that is much 
closer to the ISM value. The MMT-Pol observations re¬ 
ported here independently reconfirm the previous find¬ 
ing that the Clump is optically thick in the 1 — 3 fi m 
range (see next subsection), reaffirming the 5 x 10 -5 M 0 
minimum mass for the SW Clump and underscoring the 
need for further investigation of this peculiar hypergiant. 


3.3. Fractional Polarization versus Scattering Optical 

Depth 


Our polarimetry of IRC +10420 and VY CMa show 
their nebulae exhibit fairly high intrinsic fractional po¬ 
larizations in the infrared, with values for the various fea¬ 
tures ranging from 30% to 60% (see Table 2 summary). 
High fractional polarization is most easily reconciled with 
scattering by optically thin dust close to the plane of the 
sky. Here we assess whether the polarimetry is consistent 
with optically thin scattering. We conservatively define 
the cut-off for the transition to optical thickness at r| c > 
0.1. For reflected light from optically thin dust, the scat¬ 
tering optical depth r| c of each feature may be computed 
with: 


4t r • I% eb • 4> 

~ Ff • sin 2 (0) • <f>(0) 


(Sellgren et al. 1992). In this expression <fi is angular 
radius on the sky, I is the flux of the star at A, 0 is the 
scattering angle and 4>(0) is a phase function between 
0 and 1 that accounts for the variation with scattering 
angle 0 of the intensity of scattered light from a sphere 
with assumed optical properties. Since moderately high 
intrinsic fractional polarizations indicate scattering close 
to the plane of the sky, we adopt a scattering angle of 
0 = 90° (plane of the sky) for simplicity. Assuming 
astronomic al silicate s ph eres of size a < 0.3 /im, <f>(0 = 
90°) « 0.4 (Sh enoy et afl 2013)1. 


For each feature, we use its raw polarized intensity for 
I% eb because this is a reliable lower-limit on the total 
scattered light intensity of the feature. The resulting 
r| c lower limits discussed here are summarized in Table 
2 . For IRC + 10420’s nebula at 2.2 /im, at the peak of 
the azimuthal average of the polarized intensity at ra¬ 
dius 1.7" we find r| c ~ 0.4. For VY CMa, although the 
star is saturated in the J' images we can estimate its flux 
Fj, accurately enough in order to estimate the scatter¬ 
ing optical depth of each of the five features discussed in 
the previous subsection. We fit a Moffat profile to the 
unsaturated portion of the star’s light profile at a po¬ 
sition angle +148° E of N, which avoids faint artifacts. 
The flux Fj, obtained from this fit agrees to within ^ 
35% with the star’s flux estimated from comparison of 
the narr ow aperture photom etry (0.4" diameter) of VY 
CMa by Smith et al. (2001) with the unpolarized stan¬ 
dard star HD 22446Y. We find rj? values of 0.7 for the 
NW Arc, 1.0 for Arc 2 , 0.8 for the S Knot, and 0.5 for 
the S Arc. For the SW Clump, we find rj? > 1 and 73 % 
= 0.3. We emphasize that these are lower limits , since 
the polarized intensity of each feature is a lower limit on 
its total scattered light intensity. 

In all cases the computed lower limit to the scatter¬ 
ing optical depth is above a purely optically thin regime. 
For grains with very high albedo, it would be expected 
that multiple scatterings of the light would depolarize it 
and thus prevent the relatively high intrinsic fractional 
polarizations seen in both hypergiants’ nebulae. It has 
been shown in the visual, however, that for reflection 
from optically thick slabs relatively distant from the illu¬ 
minating star, the depolarizing effect of multiple scatter¬ 
ings on linear polari zation is less than expected at first 


■ Pc 

glance (|White||l979). The maximum linear polarization 
of scattered optical light from an optically thick slab is 
on average within a factor of 0.7 of the polarization for 
singly scattered light. White demonstrated that dust 
grains with albedos of uo < 0.4 could produce fractional 
polarizations of 40% in visual light scattered off of an 
optically thick slab. For spherical dust grains of radiu s 
a « 0.1 p m composed of astronomical silicate (Draine & 


Lee|1984 ), the albedo is about 38% at 2.2 /im. For spher- 


ical dust grains used in modeling dusty oxygen rich, late 
type stars (Suh 1999), albedos are well below 40% at 
near IR wavelengths. Our polarization observations are 
thus consistent with optically thick nebulosity for typical 
dust. 


3.4. IRC +10420: Comparison of K' Polarimetry with 
3—5 /am Images 

Here we use the image of IRC + 10420’s nebula as re¬ 
vealed by the K' (2.2 /im) polarimetry to help interpret 
the extended emission seen in the LMIRCam 3 — 5 /im 
images. We assess whether the nebula’s emission at these 
longer wavelengths is primarily scattered or thermal light 
or a combination of both. The K' polarimetry overlaps 
with the LMIRCam images starting at a radius > 0.5". 
The 3 -cr levels in the PAH 1 , PAH 2 , L' and M images lie at 
radii of 1.0", 1.2", 1.8", and 1.0" respectively. To obtain 
the broadest wavelength coverage we select an outer ra¬ 
dius of 1 . 0 " for the overlap-of-coverage region (hereafter 
OCR) that we examine. The PSF FWHM is < 0.15" 
in all four LMIRCam filters and therefore in the OCR 
the star’s light profile is negligible. This is in contrast to 
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TABLE 2 

Nebular Features’ Intrinsic Polarization and Minimum Scattering Optical Depths 


Feature 

PA+ 

(° E of N) 

distance from a 
(") 

A 

(pm) 

raw 

P A 

intrinsic p™ eb 

minimum^ r| c 

IRC +10420 

Azimuthal Avg. 


1.7 

2.2 

15% 

30% 

0.4 

VY CMa 

NW Arc 

-80 

3.4 

J' (1.3) 

35% 

35% 

0.7 

Arc 2 

-175 

4.1 

1 1 

45% 

45% 

1.0 

S Knot 

+180 

1.3 

1 1 

40% 

60% 

0.8 

S Arc 

+180 

2.4 

1 1 

40% 

60% 

0.5 

SW Clump 

-135 

1.6 

1 1 

35% 

40% 

1 

11 

-135 

1.6 

3.1 

20% 

40% 

0.3 


t Location of feature with respect to the star. 

$ Scattering optical depths are lower-limits, computed using polarized intensity as a reliable lower-limit on total scattered light intensity 
(see §3.3). 


the MMT-Pol observations at K 7 , where it was necessary 
to heavily saturate the star image in order to bring up 
the much fainter nebulosity. The OCR is depicted with 
annuli on the LMIRCam images in the bottom row of 
Figure |2] and on the K 7 image in Figure ljc). 

The nebula’s emission from K 7 througnM in the OCR 
may be a combination of thermal emission and scattered 
light. The thermal component is assumed to be unpolar¬ 
ized, while the scattered light can have both unpolarized 
and polarized components: 


jneb jneb,therm 

- 1 


jneb,scat 


A 

ir=ii 


■ f jneb,scat 

+ Oa,c/ 


1 jneb,scat \ 
+ 1 X,P ) 


(6) 

( 7 ) 


Subscript P refers to the polarized intensity, and U refers 
to the unpolarized intensity in the scattered light. We 
first consider the minimum and total scattered light flux 
at K 7 from the OCR. As in §3.1, the raw polarized inten¬ 
sity at K' (Fig. 0c)) must be scattered light from the 
nebula: I^ W P —> I^ f b p cat . This provides a lower limit 

on i 7 ^^ scat ? since the nebula must be less than 100% po¬ 
larized. Assuming that thermal emission at K 7 is negli¬ 
gible in comparison to the scattered light intensity, then 
the total scattered light intensity at K 7 is simply this 
lower limit divided by the nebula’s intrinsic fractional 
polarization: I^ scat = I™w p / p neb Per § 31) p neb „ 

.3, and therefore l k , ^ 3 x 1 K , p .the lower- 
limit and estimated total scattered intensities integrated 
over the OCR are plotted as fluxes F^, w p and Fp eb,scat 
respectively on the star’s spectral energy distribution in 
Figure [8] 

We next consider an upper limit on the contribution of 
unpolarized thermal emission to the nebular flux at K 7 . 
Given the polarimetry o f VY CMa de scrib ed in §3.2 and 
at visual wavelengths in Jones et al.| ([2007) (reproduced 
in Fig. 4b), the maximum fractional polarization of the 
scattered light from IRC -1-10420’s nebula at K 7 we would 
expect in the OCR is about 60%: 


P 


neb,scat _ 


jneb,scat 


jneb,scat 

L U 


jneb,scat — 


< 0.6 


(8) 


jneb,scat 
1 P _ 

Jneb,scat 


< 0.6 


(9) 


The nebula’s intrinsic fractional polarization of 30% 
(§3.1) is determined by the nebula’s total scattered plus 
thermal emission: 


P 


neb 


jneb,scat 

_Tp_ 

Jneb,scat _|_ Jneb,therm 


= 0.3 


( 10 ) 


Comparing the maximum and intrinsic polarizations, we 
can write: 


jneb,therm jneb,scat 


(ii) 


Thus, the maximum contribution of unpolarized thermal 
emission to the observed flux from the OCR in the K' 
filter would be roughly equal to the scattered light con¬ 
tribution. 

The intensity of scattered light at a given wavelength is 
directly proportional to the illuminating star flux, which 
we designate as and the grains’ albedo uj\. We scale 
jneb,seat f rom j£' to the four longer wavelengths A = 
PAH1, PAH2, L' and M using ratios for the illuminat¬ 
ing flux and grain albedo at the two wavelengths. We 
are assuming the OCR is outside of the region produc¬ 
ing the bulk of the emission that makes up the SED of 
IRC+10420. In this way we can take the observed SED 
as the illuminating flux F£ and express the following re¬ 
lationship: 


jneb,scat jneb,scat 

7 A J K 


F* 

r K 


+ ) ( ^ ) (12) 


Note that the effects of interstellar reddening on the 
star’s SED and the nebula are the same, so the extrap¬ 
olation from K 7 to the longer wavelengths can be done 
directly with the reddened (observed) SED. The ratio of 
the observed fluxes is F^/F^, ~ 2 for all four LMIRCam 
filters. 

Assuming silicate grains of radius a < A (e.g., a < 0.1 
H m) the ratio of the albedos lo\/lok’ for the four filters 
are 0.39, 0.36, 0.26, and 0.13. These albedo ratios are 
typical ra tios for astro nomical silicates for the interstellar 
medium (Dr aine fc Lee|1984) and oxygen rich mass-loss 
winds (fSuhf|1999|). With these values in Equation (12), 
we predict the 3 — 5 fi m scattered light fluxes depicted 
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with open blue circles in Figure [8j On the SED in Figure 
[8j the nebula’s integrated 3 — 5 /im fluxes for the OCR 
Rom 0.5" to 1.0 " are plotted as grey circles. The open 
red circles are the excess flux in each of these 3 filters 
(total — predicted scattered). As can be seen in Figure 
[8j the total nebular emission at 3 — 5 /im is about a fac¬ 
tor of 10 brighter than scattered light alone can account 
for. It is most likely the nebula’s emission at ^ 3 /im 
and beyond is primarily thermal in origin. However, we 
have just shown that at most half of the flux at 2.2 /im 
can be thermal, and this puts strong constraints on any 
explanation for the surprisingly large 3 — 5 /im fluxes. 

On the SED in Figure [8] we plot the spectrum of a 500 
K emissivity-modified blackbody which passes through 
the L' excess flux. This choice of temperature was cho¬ 
sen to prevent too much thermal flux at 2.2 /im (a higher 
temperature) and too much flux in the mid-IR (a lower 
temperature) which would compete with or exceed the 
observed mid-IR flux from the entire system. A lower 
temperature such as 400 K, for example, would require 
fully half of the entire 10 /im silicate feature to arise 
from the OCR alone, which seems unlikely. The black- 
body equilibrium temperature is 120 K, assuming L^ ol = 
5 x 10 5 Lq and Dm 5 kpc. Factoring in wavelength- 
dependent emissivity Q^ bs to account for less efficient 
emission in the infrared, the dust equilibrium tempera¬ 


ture is T eq = ((Quv,vis) / {Qir}) 1 ^ 4 T bb , where (Quv,vis) 
and (Qir) are the Planck-averaged emissivities for the 
wavelength ranges where the grains absorb and emit, re¬ 
spect ively^^oi_thety2iva]_^ilicate grains we have been 
using (Draine & Lee|1984 Suh 1999) this factor raises the 
equilibrium temperature to ~ 210 K. Although warmer 
than blackbody equilibrium, this is still too cool to be 
compatible with the observed mid-IR flux. 

The ^ 500 K color temperature we require is 2 to 
3 times higher than the expected equilibrium tempera¬ 
ture. One possibility is emission from transiently heated 
small dust grains, which can temporarily reach much 
higher equilibrium temperatures from absorption of a 
single stellar photon with an energy comparable to the 
grain’s heat capacity. These grains are usually associ¬ 
ated with reflection ne bulae illuminated by a source of 
near-UV photons, e.g., Sellgren et al. (1983). However, 
the spectrum of emission from these very small grains 
corresponds to color temperatures much higher than the 
500 K we requi re to ex plain the emission from the OCR. 
For example, Draine & Li (2007) compute SEDs for a = 
5 A silicate grains exposed to the interstellar radiation 
field (roughly equivalent to an F0 star, similar to the 
IRC+10420). They find such a tiny grain’s spectrum to 
be relatively flat from 2.2 — 5 /im (see their Figure 10), 
which is incompatible with our observations. We do note 
that our 500 K color-temperature is close to with the gas 
temperatur e at this radius (r ~ 0.75" = 5.6 x 10 16 cm) 
modeled by Castro-Carrizo et al. (2007) in their study of 
CO emission around IRC +10320. These authors mod¬ 
eled the CO emission as arising from two spherical shells, 
with the inner shell spanning from 0.3" to 1.7" in radius, 
very similar to our OCR. At the 0.75" mean radius of 
the OCR we examine, their best-fit model predicts a gas 


temperature of T eq = 460 K. While our polarimetry ob¬ 
servations clearly require a nebular material lying largely 
in the plane of the sky as opposed to a shell, the similarity 
in their gas temperature and the dust color temperature 
we derive indicates that higher than equilibrium temper¬ 
atures are physically real. 


4. CONCLUSIONS 

1. We present images of the nebulosity surrounding the 
hypergiant VY CMa in polarized light with a resolution 
of 0.2" at 1.3, and 3.1 /im. Many nebular features are 
seen in polarized intensity at 1.3 fi m and one feature, the 
SW Clump, is very prominent at 3.1 /im as well. 

2. Extended emission around VY CMa shows both 
high fractional polarization and high scattering optical 
depths. The surface brightness is consistent with scat¬ 
tered light alone. The high fractional polarization is con¬ 
sistent with grain albedos considerably less than 1. The 
required albedos are easily compatible with the typical 
silicate grains used to model dusty stars. 

3. The polarimetry of VY CMa independently con¬ 
firms that its Southwest (SW) Clump is optically thick 
from 1.3 — 3.1 /im. This affirm s the minim um mass of 
the SW Clump computed by Shenoy et al. (2013) from 
its scattered light in the infrared. Absence of thermal 
dust emission at mm wavelengths from the SW Clump 
presents a puzzle. 

4. We present images of the nebulosity surrounding 
IRC+10420 in polarized light at 2.2 /im with a resolution 
of 0.2" and total flux at 3.3, 3.4, 3.8, and 4.9 /im with a 
resolution of 0.1" — 0.15". 

5. The polarized intensity image of IRC+10420 at 
2.2 /im shows a relatively uniform nebula largely in the 
plane of the sky extending from 0.5" — 2.5" radius (2500 
— 12500 AU for D = 5 kpc). The surface brightness 
of this low-latitude ejecta is compatible with optically 
thick scattering and, similar to VY CMa, grain albedos 
compatible with typical astronomical silicates. 

6. In the 3 — 5 /im band, images of IRC+10420 show a 
strong extended component that overlaps with the neb¬ 
ula seen in polarized intensity at 2.2 /im. The flux 
from this extended component is an order of magnitude 
brighter than can be explained by simple extrapolation of 
the scattered light seen at 2.2 /im. We hypothesize grains 
warmed to a temperature higher than the expected grain 
equilibrium temperature, but consistent with the local 
gas temperature in this region. 
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Fig. 8.— Spect ral energy distribution of IRC +10420. The black dots are point-source photometry compiled from Jones et al. (1993}, 
|Oudmaijer et al. (1996} and the IRAS and MSX Point Source Catalogs (with color correction). The thin black line is the ISO SWS 
spectrum. is the polarized intensity at A = K 7 (2.2 p,m) integrated over the overlap-of-coverage region (OCR, the annulus depicted 

in Figure [jpc) and the bottom row of Figure [5}. This is a lower limit on the scattered light flux at K' from this region. The filled blue 
circle F™ e , ,scat is the estimated total (unpolarized + polarized) scattered light flux for ~ 0.3 (§4.1). The blue dashed line scales this 
total scattered light flux to A = PAH1 (3.29 /i m), PAH2 (3.40 jt/m, L' (3.83 n m) and M (4.9 /i m) using Equation (12) in order to predict 
the nebula’s scattered light flux at those wavelengths (open blue circles). The filled grey circles are the nebula’s total 3 — 5 fi m fluxes in 
those filters, with uncertainty smaller than the symbol size for all except the PAH1 filter. The open red circles are the excess flux (total — 
scattered). The red dashed line is an emissivity-modified A B\(T eq ) curve for T eq = 500 K. 
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